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Abstract 

; , Musculoskeletal disorder (MSD) is one of the major health problems in mechanical work especially in manual handling jobs. 
^ Muscle fatigue is believed to be the main reason for MSD. Posture analysis techniques have been used to expose MSD risks of the 
I I work, but most of the conventional methods are only suitable for static posture analysis. Meanwhile the subjective influences from 
the inspectors can result differences in the risk assessment. Another disadvantage is that the evaluation has to be taken place in the 
workshop, so it is impossible to avoid some design defects before data collection in the fleld environment and it is time consuming. 
In order to enhance the efficiency of ergonomic MSD risk evaluation and avoid subjective influences, we develop a new muscle 
fatigue model and a new fatigue index to evaluate the human muscle fatigue during manual handling jobs in this paper. Our new 
fatigue model is closely related to the muscle load during working procedure so that it can be used to evaluate the dynamic working 
process. This muscle fatigue model is mathematically validated and it is to be further experimental validated and integrated into 
a virtual working environment to evaluate the muscle fatigue and predict the MSD risks quickly and objectively. 
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Relevance to industry 

On : 

' Muscle fatigue is one of the main reasons causing MSDs in industry, especially for mechanical work. Correct evaluation of muscle 
^ fatigue is necessary to determine work- rest regimens and reduce the risks of MSD. 
• ^ 

$^ Key words: Muscle fatigue; muscle fatigue index; muscle fatigue model; virtual reality; virtual environment 



1. Introduction 

Musculoskeletal disorder is defined as injuries and dis- 
orders to muscles, nerves, tendons, ligaments, joints, 
cartilage and spinal discs and it does not include in- 
j uries resulting from sl i ps, tri ps, falls or similar accidents 
("Maier and Ross-Motal. bood ). From the report of HSE 
( HSE , 2005) and 199 3 WA State Fund comp e nsable claims 

over 



200,^ . 



(|State of Washington Department of Laboii . 
50% of workers in industry have suffered from muscu- 
loskeletal disorders. According to the analysis in the book 
Occu pational Biomechanics ( Chaffin and Andrerssoil 
overexertion of muscle force or frequent high mus- 
cle load is the main reason to cause muscle fatigue, and 
further more, it results in acute muscle fatigue, pain in 
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muscles and the worst functional disability in muscles and 
other tissues of human body. Hence, it becomes an impor- 
tant mission for ergonomists to find an efficient method to 
evaluate the muscle fatigue and further more to decrease 
MSD risks caused by muscle fatigue. 

There are various tools available for ergonomists to eval- 
uate the MSD risks, most of th em are listed and compared 
in paper (jLi and Bucklel . ll999[l . For general posture analy- 
sis, Posturegram, Ovako Working Posture Analyzing Sys- 
tem (OWAS), posture targeting and Quick Exposure Check 
(QEC) were developed. Further more, some special tools 
are designed for specified parts of human body. For exam- 
ple, Rapid Upper Limb Assessment (RULA) is designed for 
assessing the severity of postural loading and is particularly 
applicable for sedentary job s. The similar system s include 



HAMA, PLIBEL and so on (jStanton et al.l . i2004h 



In these methods, several "risk factors" of mechanical 
jobs are taken into consideration, like physical work load 



factors, psychosocial stressors and individual factors. The 
level of exposure to physical work load can be evaluated 
with respe ct to intensity (or mag nitude) , repetitiveness and 



duration (|Li and Buckld . Il999f l. Using these tools, MSD 



risks can be effectively reduced, but there are still several 
limitations. At first, even just for lifting job, five prevailing 
tools (NIOSH Hfting index, ACGIH TLV, 3DSSPP, WA L 
& I, Snook) were developed. The evaluation results of them 
for a sa me task were different, and sometimes even contra- 
dictory Russell et al.l . 120071 ) . Second, most of the tradi- 
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tional methods have to be carried out in field environment, 
there is no immediate result from the observation and it is 
time consuming for the repetition. Further more, there are 
different evaluation results from different subjects even us- 
ing a same method ( Lamkull et al. . 20071 ). At last, only pos- 
ture information and limited working conditions are con- 
sidered in these methods, and thus they are not able to 
estimate the MSD risks into detailed analysis and can not 
further enhance the safety of the work. The most significant 
problem is that muscle fatigue prediction and assessment 
tool is still in blank in conventional methods. 

In order to evaluate the human work objectively and 
quickly, virtual human techniques have been developed 
to facilitate t h e ergq nomic evaluati on. For examp l e. Jac k 
(iBadler et all. Il993h. E rgoMan (ISchaub et al.L 
3DSSPP ("ChafRd ll969|), Santos 

IVignes . ■2004f ) and so on, these human modeling tools 
are used in field of automotive, military, aerospace and 
industrial engineering. These human modeling tools are 
mainly used for visualization to provide information about 
body pos ture, reachabili t y and field of view and so on. 
In paper (^Javar am et al. . 20061 ) . a method to link virtual 
environment (Jack) and quantitative ergonomic analysis 
tools (RULA) in real time for occupational ergonomics is 
presented, and it allows that ergonomic evaluation can be 
carried out in real time in their prototype system. But 
until today, these is still no fatigue index available in these 
virtual human tools for dynamic working process. It is 
necessary to develop the muscle fatigue model and then 
integrate it into the virtual human softwares to evaluate 
the muscle fatigue and analyze the mechanical work into 
details. 

For objectively predicting muscle fatigue, several mus- 
cle fatigue models and fatigue index ha ve been pro posed in 
publications. In a s e ries of p u blicati ons (jWexler et al. . 1997: 
Ding et al.l . [2000al lbl. l2002bl lal l2003f) . Wexler and his colleges 
have proposed a new muscle fatigue model based on Ca^+ 
cross-bridge mechanism and verified the model with stimu- 
lation experiments, but it is mainly based on physiological 
mechanism and it is too complex for ergonomic applica- 
tion. Further more, there are only parameters available for 
quadriceps. This model was integrated into VSR system for 
several muscle fatigue simulations w ith lots of lim itations 
jVignes, 2004). Taku Komura et al. (jKomura et al.. 1999, 
2000f) have empl oyed a muscle fatig ue model based on force- 
pH relationship f Ciat et al , 1993f ) in computer graphics to 
visualize the muscle fatigue, but it just evaluated the mus- 



cle fatigue at a time instant and cannot evaluate the overall 
muscle fatigue of the working processes. Meanwhile in this 
pH muscle fatigue model, although the force generation ca- 
pacity can be mathematically analyzed, the influences from 
the muscle forces are not conside red enough. Rodriguez 
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propos es a half-joint fatigue model (jRodrfguez et al 
2003allbl ). more exactly a fatigue index, based on mechan- 



ical properties of muscle groups. The fatigue is quantified 
by accumulation of the proportion of actual holding time 
over the predicted maximum holding time. With this half- 
joint model, it can adjust human posture during a working 
process dynamically, but it cannot predict individual mus- 
cle fatigue due to its half-joint principle. The maximum 
endurance equation of this model was from static posture 
analysis and it is mainly s uitable for evaluating static pos- 



tures. In (|Liu et al.l . l2002f ). a dynamic muscle model is pro- 



posed based on motor units principles, but there are just 
parameters available under maximum voluntary contrac- 
tion situation which is rare in manual handling work. 

In this paper, we are going to propose a dynamic mus- 
cle fatigue model and a fatigue index with consideration 
of muscle load history and personal factors. This fatigue 
model is going to be verified with comparison of the pre- 
vious static endurance time models and several dynamic 
muscle fatigue models. At last, we are going to propose the 
experimental validation procedure in a virtual environment 
framework. 



2. Dynamic Muscle Fatigue Model 

We believe that the fatigue of muscle is closely related 
to the external load of the muscle with the time and the 
strength of the muscle. These factors can represent the 
physical risk factors mentioned before: the external load the 
muscle with the time can include the intensity (or magni- 
tude) , repetitiveness and duration; and the muscle strength 
can be determined individually. Thus, the muscle force his- 
tory and maximum voluntary contraction (MVC) are taken 
to construct our muscle fatigue model. MVC is defined 
as "the force generated with feedback and encouragement, 
when the subject believes it is a maximal effort" (VQllcstadl, 
Il997h . The effect of maximum voluntary contraction on en- 
durance time is often used in erg onomic appli cations to 
define the worker capabilities (jOarg et al.l . 12002 ). Here, we 
are going to use MVC to describe the maximum force gen- 
eration capacity of an individual muscle. 

Our fatigue index is trying to describe the human feelings 
(subjective evaluation) about the fatigue and it is expected 
to have a close correlation between subjective evaluation 
and objective evaluation about muscle fatigue during me- 
chanical work. 

For ergonomic evaluation, there are several self-report 
methods to assess the physical load, body discomfort or 
work stress. These subjective assessments of body strain 
and discomfort have been the most frequently used form 
due to the ease of use and apparent face validity. But subjec- 



tive ratings are prone to many influences. This kind of ap- 
proach has too low vahdity and rehability. But ergonomists 
have to concentrate themselves on the feeling of the work- 
ers. Several authors even insist that "If the person tells you 
that he is loaded and effortful, then he is loaded and ef- 
fortful whate ver the behavioral and performances measures 
may show" (|Li and Buckl i,[l99i). Thus, to combine the 
subjective evaluation and objective evaluation can reduce 
the contradictory problems between them and can rate the 
muscle fatigue correctly. 

The general feeling of human body about fatigue is: the 
larger the force is, the faster people can feel fatigue; the 
longer the force maintains, the more fatigue; the smaller 
capacity of the muscle is, the more easily we can feel fa- 
tigue. Based on this description, the equation [T] is used to 
describe the subjective evaluation. The parameters used in 
the equations are listed and described in table[TJ This equa- 
tion can be explained as follows: 

(i) Fcem describes the capacity of the muscle after some 
contraction, the rest force generation ability. 

(ii) Fioad{t) / Fcem{t) is relative load which describes the 
muscle force relative to the capacity of the muscle at 
a time instant t. 

(iii) MVC / Fcemit) describes the smaller capacity of the 
muscle, the faster the muscle get fatigued. 

Table 1 

Parameters in Dynamic Fatigue Model 



Item 


Unit 


Description 


MVC 


N 


Maximum voluntary contraction, maximum ca- 
pacity of muscle 


Fcem {^) 


N 


Current exertable maximum force, current capac- 
ity of muscle 


Fload{t) 


N 


External load of muscle, the force which the mus- 
cle needs to generate 


k 


•min~^ 


Constant value, 1 


U 


min 


Fatigue Index 


%MVC 




Percentage of the voluntary maximum contraction 


Imvc 




%MVC/100 



dU{t) MVC Fioadit) 



(1) 



Meanwhile, the current maximum exertable force Fcem 
is changing with the time due to the external muscle load. 
It makes sense that: The larger the external load, the faster 
Fcem decreases; the smaller Fcem is, the more slowly Fcem 
decreases. The differential equation for Fcem is equation[2l 



dFcemit) jFcemit) , . 

The integration result of equation [5] is equation [31 



(2) 



P Fload{u) 
I — k du 

Fcemit) = MVC eo ^^^^ (3) 



MVC 



du 



(4) 



MVC is a constant value for an individual person, so we 
can change the equation[3]into equation^ If Fioad/MVC is 
constant and equals to C, then F{t) = Ct, equation [3] can 
be further simplified. This constant case can occur during 
static posture and static load. 



Fcem jt) _ up(^i^ 



-kCt 



(5) 



MVC 

The feeling of our fatigue is a function below and which 
is closed related to M VC and Finadit). MVC can rep resent 
the personal factors ( Chaffin and Andrersson . 1993), and 



Fioad{t) is the force exerted on the muscle along the time 
and it reflects the influences of external loads equation [6l 



U{t) = 



1 



,2kF{t) 



1 



,2feF(0) 



(6) 



2k 2k 
In this model, personal factors and external load history 
are considered to evaluate the muscle fatigue. It can be 
easily used and integrated into simulation software for real 
time evaluation especially for dynamic working processes. 
This model still needs to be mathematically validated and 
ergonomic experimental validated, and meanwhile muscle 
recovery procedure should also be included to make the 
fatigue index completed. 

3. Static Validation 

3.1. Validation result 

Our dynamic muscle fatigue model is simply based on 
the hypothesis on the reduction of the maximum exertable 
force of muscle. It should be able to describe the most spe- 
cial case - static situations. In static posture analysis, there 
is no model to describe the reduction of the muscle capacity 
related to muscle force, but there are several models about 
maximum endurance time (MET) which is a measurement 
related to static muscular work. MET represents the max- 
imum time during whi ch a static load can be maintained 



jElahrache et al.ll2006l ). The MET is most often calculated 



in relation to the percentage of the voluntary maximum 
contraction {%MVC) or to the relative force {fhivc = 
%MVC /lO O) required by th e task. These models which are 
cited from (jElahrache et al.i . 2006) are listed in table [2l 

In our dynamic model, suppose that Fioad{t) is constant, 
and it represents the static situation. MET is the duration 
in which Fcem falls down until the current Fioad- Thus, 
MET can be figured out in equation [7] and [H 



,{t) = MVCeo 



f-k- 



o.d{u) 



MVC 



du 



= Floadit) (7) 



Assume that F{t) is: 



t = MET = -- 



, Fioad{t) 
Pload{t) 



k Jmvc 



Intraclass Correlation between Dynamic model and General models 



(8) 



MVC 

In order to analyze the relationship between MET of 
our dynamic model and the other models, two correlation 
coefficients are calculated. One is Pearson's correlation r in 
equation [HI and the other one is intraclass correlation ICC 
in equation llOl r indicates the linear relationship between 
two random variables and ICC can represent the similarity 
between two random variables. The closer r is to 1, the more 
the two models are linear related. The closer ICC is to 1, 
the more similar the models are. M Si 



betu 



is the mean 

square between different MET values in different Jmvc 
values, MSwithin is the mean square within MET values in 
different models at the same Jmvc level, p is the number 
of models in the comparison. In our case, we compare the 
other models with our dynamic model one by one, k equals 
to 2. The calculation results are shown in table[2]and figure 

[BtolHl 



ICC = 



ZiA, ^ A){Br, - B) 

n 

E(A„-A)2^(B„-B)2 

71 

MSbet ween 

Endurance Time in General Models 



(9) 



(10) 
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Fig. 1. Endurance time in general models 



3.2. Discussion 

From the comparison results of the static validation, it is 
obvious that MET model derived from our dynamic model 
has a great linear correlation with the other experimental 
static endurance models, and almost all the Pearson's cor- 
relation r are above 0.97. 

From ICC column, the similarity between our dynamic 
model and the other models in descend sequence is: elbow 
models, hand model, general model, shoulder models and 




2 3 4 5 6 7 

Endurance Time in Dynamic model [min] 



Fig. 2. ICC of general models 



Endurance Time in Shoulder Models 




Fig. 3. Endurance time in shoulder endurance models 



intraclass Correlation between Dynamic model and Shoulder models 




2 3 4 5 6 7 

Endurance Time in Dynamic model [min] 



Fig. 4. ICC of shoulder endurance models 

hip/back model. For elbow models, the average ICC is ap- 
proximate 0.9, but for the back/hip models, ICC varies 
from -0.057 to 0.9447. The explanation is: in shoulder and 
back/hip of human body, the anatomical structure is in a 
much more complex way than in the elbows and hands. 
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Table 2 

Static validation results telahrache et all [20o3) 



Model 


Equation 


r 


ICC 


General models 



Rohmort MET ■- 

Monod and Scherrer MET ■- 

Huijgens MET 

Sato et al. MET -- 

Manenica MET ■- 

Sjogaard MET -- 

Rose et al. MET -- 



-1.5 + 



fMVC fl,^c fl 

0.4167 (/mvc -0.14)-2-4 

0.3802 (/mvc -0.04)-i-44 

14.88 exp(-4.48/MVc) 
0.2991 /-//p* 

7.96 exp(-4.16/MVc) 



0.9937 0.8820 

0.8529 0.6474 

0.9964 0.8800 

0.9992 0.8512 

0.9927 0.9796 

0.9935 0.9917 

0.9897 0.7080 



Upper limbs models 



Shoulder 
Sato et al. 
Rohmert et al 



' MVC 
-1.658 



MET = 0.398/^ 
MET = 0.2955 f-^y^ 
Mathiassen and Ahsberg MET = 40.6092 exp(-9.7/MVC 



Cars 



MET = 0.5618/" 



0.9997 0.7188 

0.9987 0.5626 

0.9783 0.7737 

0.9981 0.9029 



Elbow 



Hagberg 


MET = 


0-298 /l/v-c 


0.9935 


0.9921 


Manenica 


MET = 


20.6972 exp(-4.5/MVc) 


0.9929 


0.9271 


Sato et al. 


MET = 




0.9838 


0.9712 


Rohmert et al. 


MET = 


0.2285 


0.9997 


0.7189 


Rose et al.2000 


MET = 


20.6 exp(-6.04/Mvc) 


0.9986 


0.9594 


Rose et al.l992 


MET = 


10.23 exp(-4.69/MVc) 


0.9943 


0.7843 


Hand 










Manenica 


MET = 


20.6972 cxp(-4.5/A/vc) 


0.9929 


0.9840 


Back/hip models 


Manenica (body pull) 


MET = 


27.6604 exp(-4.2/MVc) 


0.9901 


0.6585 


Manenica (body torque) 


MET = 


12.4286 exp(-4.3/MVc) 


0.9911 


0.9447 


Manenica (back muscles) MET = 


32.7859 exp(-4.9/Mvc) 


0.9957 


0.7306 


Rohmert (posture 3) 


MET = 


0.3001 


0.9745 


0.5353 


Rohmert (posture 4) 


MET = 


1.2301 


0.9989 


0.7041 


Rohmert (posture 5) 


MET = 


3.2613 /"V'^s 


0.9984 


-0.057 



In this case, during these experimental models, the mea- 
surement of MVC is an overall performance of the muscle 
group, but not an individual muscle. Meanwhile from fig- 
ure [71 the differences between the experimental models for 
hip/back are much greater than in the other models like in 
elbow models (figure^. It can be explained as: in different 
working conditions (for example, different postures), the 
engagement of the muscles in the task in the hip/back of 
human body varies due to the complexity of the structure. 
Thus, the error between our model and the other models 
is closed related to the complexity of the structure and our 
model can fit most of the experimental models with a high 
similarity. 

Overall, the static validation can prove that our dynamic 
model can be used to predict MET in static situations. 



4. Dynamic Validation 

4.1. Validation result 

Static validation results have shown that our dynamic 
model can be used to predict the MET for general static 
load and even for some specific body parts. But static proce- 
dures are still quite different from dynamic situations, thus 
our dynamic model need to be under examination with the 
other dynamic models. For this objective, we are going to 
verify our dynamic model through comparison with some 
existing muscle fatigue models, quantitatively or qualita- 
tively. 

In paper (jFreund and Takalal . l2002l) . a muscle fatigue 



Endurance Time in Elbow Models 



Intraclass Correlation between Dynamic model and Hip/back models 
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Endurance Time in Dynamic model [min] 



Fig. 5. Endurance time in elbow endurance models 



Fig. 8. ICC of hip/back models 



Intraclass Correlation between Dynamic model and elbow models 



1 = 





Hag berg 




Manenica 




Sato 




Rohmert 


+ 


RoseOO 




Rose92 




Dynamic 
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Endurance Time in Dynamic model [min] 



Fiff. 6. ICC of elbow endurance models 



Endurance Time in Hip/Back Models 



Manenica pull 

Manenica torque 

Manenica muscle 

Rohmert3 

Rohmert4 

Rohmert5 
* - Dynamic 
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Fig. 7. Endurance time in hip and back models 

model was proposed and integrated into a dynamic model 
of forearm. In this model, the muscle was treated like a 
kind of reservoir, and force production capacity reduces 
with the time that the muscle is contracted. varies be- 
tween and the upper limits of the muscle force 5' . In this 



model, the recovery and decay rates depend on S*' — and 
muscle force S (equation [Tl]). The constants a and [3 were 
obtained by fitting the solution using experimental results 
from static endurance time test. In this model, muscle force 
is taken into consideration as a factor causing muscle fa- 
tigue, and further more, muscle force production capacity 
S*" was proposed just like in our dynamic model fcem to 
describe the capacity of the muscle after performing cer- 
tain task. But in this model, the force production capacity 
and the muscle force are decoupled with each other which 
is different in our model. 



dt 



(11) 



Wexler's dynamic muscle fatigue model based on Ca?'^ 
cross-bridge mechanism can also verify our dynamic model 
quantitatively. This model can be used to predict the mus- 
cle force fatigue under different stimulation frequencies. 
From figure [H it is clear that the faster the stimulation 
frequency is, the larger the force can be generated by the 
muscle. If each curve is normalized into figure [TOl it is obvi- 
ous that: the larger the peak force is, the faster the muscle 
gets fatigue. But this model can only be used to predict the 
muscle force under a stimulation pattern and it can not pre- 
dict or evaluate the muscle force generation capacity under 
a stimulation pa ttern when the m uscle is getting fatigue. 

In the paper (jLiu et al. . l2002h . another dynamic model 
of muscle activation, fatigue and recovery was presented. 
This model is based on biophysical mechanisms: a muscle 
consists of many motor units which can generate force or 
movement. The number of the motor units depends on the 
size and function of the muscle. The generated force is pro- 
portional to the activated motor units in the muscle. The 
brain effort B, fatigue property F and recovery property R 
of muscle can decide the number of activated motor units. 
The relationship is expressed by equation II 21 The param- 
eters in this equation is explained in table [3l 



Frequency 100 [ms] 
Frequency 50 [ms] 
Frequency 10 [ms] 
Frequency 2 [ms] 




Time [ms] 



Fig. 9. Maximum exertable force and time relationship in Wexler's 
Model 



Frequency 100 [ms] 
equency 50 [ms] 
equency 10 [ms] 
equency 2 [ms] 




Time [ms] 



Fig. 10. Normalized maximum exertable force and time relationship 
in Wexler's Model 



M 



M. 
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Fig. 11. Illustration of the active motor model ihiu et al.l |200 



dMA 
dt 

Muc 



= BMuc-FMa^ 
= FMa ~ RMf 
= Mo -Ma- Mf 



RMf 



(12) 



When t — Q under the initial conditions of Ma — 0, 
Mf — 0, Muc — Mq, we can have equation [T3l 



MAjt) 
Mn 



1 + 7 (l + 7)(/3-l-7) 
_ P-1 -0Ft 



-(l+7)Ft 



(13) 



Table 3 

Parameters in Active Motor Model 



Item Unit Description 



F fatigue factor, fatigue rate of motor units 

R s^^ recovery factor, recovery rate of motor units 
B s^^ brain effort, brain active rate of motor units 
total number of motor units in the muscle 
number of activated motor units in the muscle 
number of fatigued motor units in the muscle 
number of motor units still in the rest 
B/F 
R/F 



Mo 
Ma 
Mp 

Muc 

/3 
7 



In our fatigue model, we assume that there is no recovery 
during mechanical work, and the workers are trying their 
best to finish the work which means the brain effort is quite 
large. In this assumption, we set 7=0 and (3 — s- 00, then 
the equation [14] represents the motor units which are not 
fatigued in the muscle and Ma (t) /Mq represents the muscle 
force capacity. We can simplify the equation[T3]to equation 
1141 which does have the same form of our dynamic model 
equation [SI 



MAjt) 

Mn 



-Ft 



-l3Ft 



-Ft 



(14) 



This fat i gue m odel has been experimentally verified 



I Liu et al. . 2002f) . In the experiment, each subject per 



formed an MVC of the right hand by gripping a handgrip 
device for 3 min. And the fitting curve from the experiment 
result has almost the same curve of our model in MVC 
condition (figure [T2|) . In this model, F and R are assumed 
to be constant for an individual under MVC working con- 
ditions. There is no experiment result for F and R under 
the other load situations, thus this muscle fatigue model 
can only verify our model in MVC condition. 

4.2. Discussion 

Through dynamic validation, our dynamic model is 
quantitatively and qualitatively verified with the other 
existing muscle fatigue models. The fatigue model used in 
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comparison result between active motor model and dynamic model 



- Active Motor Experiment result,F=0. 0206, |3=254,y=0. 398 

- Dynamic Model, f /MVC = 1 



20 30 
time unit [s] 



Fig. 12. Comparison between the experimental result of the active 
motor model and dynamic model in the maximum effort 

the forearm used the same conception Hke in our fatigue 
model: the muscle force capacity is related to muscle force 
with time. Wexler's model based on Ca^+ cross-bridge 
shows the reduction of the muscle force during the time 
under different stimulation frequencies, the reduction of 
the muscle capacity shows the same trend like in our mus- 
cle fatigue model. With comparison of the active motor 
model, the muscle force can be expressed in the same form 
under extremity situation. But in the active motor model, 
only parameters are available for MVC contraction case. 
The active motor does not supply further validation for 
other load situations. But we believe that the muscle load 
can influent the fatigue factor in the active motor model. 

Overall, our dynamic model is simple and easy to use 
and it can evaluate the muscle fatigue during a dynamic 
working process. 

5. Experimental Validation 
5.1. Objective and methodology 

After mathematical validation of our dynamic model, we 
are going to construct a virtual reality framework to verify 
the fatigue index in experimental environment. The objec- 
tive of the experimental validation is to find the correlation 
between subjective fatigue evaluation results and objective 
fatigue evaluation results. In our hypothesis, if both results 
are highly linear related, that means the objective meth- 
ods can represent the subjective feeling of fatigue and then 
it can be further integrated into virtual human simulation 
software to evaluate human muscle fatigue during manual 
handling work. 

In order to realize our objective, two evaluation systems 
should be established: subjective evaluation system and ob- 
jective evaluation system. Subjective evaluation methods 
have been mentioned in many papers, and here we are fo- 
cusing on objective evaluation system. For objective evalu- 
ation, according to our model, muscle force at each time in- 



stant during the work should be calculated. It requires dy- 
namic information like acceleration, velocity, position and 
external load. These information can be further input into 
dynamic formulas to calculate the muscle force. In order to 
get these information, motion capture method should be 
involved. As mentioned, the traditional methods have sev- 
eral drawbacks for ergonomic analysis, a virtual environ- 
ment needs to be constructed to decrease the time cost and 
physical prototype cost. 



5.2. System Structure 

The overall objective of the framework is to evaluate hu- 
man work and predict potential human MSD risks dynam- 
ically, especially for human muscle fatigue. The function 
structure of the framework is shown in figure [T3] and dis- 
cussed below. 

In order to avoid field-dependent work evaluation, vir- 
tual reality techniques and virtual human techniques are 
used. Immersive work simulation system should be first 
constructed to provide the virtual working environment. 
Meanwhile, virtual human should be modeled and driven 
by the motion capture data to map the real working pro- 
cedure into the virtual environment. Haptic interfaces can 
be used to enable the interactions between the worker and 
virtual environment. 

For any ergonomic analysis, data collection is the first 
important step. All the necessary information needs to be 
collected for further processing. From section [2l necessary 
information for dynamic manual handling jobs evaluation 
consists of motion, forces and personal factors. Motion cap- 
ture techniques can be applied to achieve the motion infor- 
mation. In general, there are several kinds of tracking tech- 
niques available, like mechanical motion tracking, acoustic 
tracking, magnetic tracking, optical motion tracking and 
inertial motion tracking. Each tracking technique has its 
advantages and drawbacks for capturing the human mo- 
tion. Hybrid motion tracking techniques can be taken to 
compensate the disadvantages and achieve the best motion 
data. 

Force information can be recorded by haptic interfaces. 
Haptic interface is the channel via which the user can com- 
municate with virtual objects through haptic interactions, 
and the interaction data between the worker and the vir- 
tual environment are also significant for evaluating other 
ergonomic aspects. Individual factors can be achieved 
from anthropometrical database and some biomechanical 
database. 

All the information, such as motion information, force 
history and interaction events, is further processed into Ob- 
jective Work Evaluation System (OWES). The output of 
the framework is evaluation results of the mechanical work. 
There are many ergonomic aspects of a mechanical work. 
For each aspect, corresponding criteria should be estab- 
lished to assess the dynamic work process. Further more, 
these criteria can also be applied just into commercial hu- 
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man simulation software to generate much more naturally 
and realistically human simulation. 
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Fig. 13. The framework of the dynamic work evaluation system 

Three hardware systems are employed to construct a 
prototype system to realize this framework: virtual sim- 
ulation system, motion capture system and haptic inter- 
faces. Virtual Simulation system consists of graphic simu- 
lation module and display module. Simulation module exe- 
cutes on a computer graphic station, and display module is 
composed of projection system and head mounted display 
(HMD). Simulation module is in charge of graphic process- 
ing and display control. The projection system and HMD 
system can visualize the immersive environment. Motion 
Capture system is taken to capture the movement of the 
worker to collect the real working data. The Optical Motion 
Cap ture System de veloped by Tsinghua University VRLIB 
( Wa ng et al. . 20061 ) is used, and it can work at 25 Hz data 
update rate. Each key joint of human body is attached by 
an active marker, and the overall skeletal structure of hu- 
man body can be represented by 13 markers. Meanwhile, 
data gloves are used to capture the movement of fingers. 
Haptic system is also going to be integrated into this pro- 
totype, and now we are still under preparation. 
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Fig. 14. Scheme of Hardware system for the experimental vaUdation 



6. Conclusion and Prospects 

In this paper, we present a new muscle fatigue model and 
a new fatigue index based it. The mathematical validation 
has proved that the new model is very simple and easy to 
predict the MET in static postures. The dynamic validation 
result proves that the new model can be further extended 
into dynamic working process to predict the muscle fatigue. 
For experimental validation, a virtual reality framework is 
under construction to verify the linear relation between our 
fatigue index and subjective evaluation result. 

The prototype system cannot only evaluate the fatigue, 
but also the other aspects of work by extending the evalua- 
tion criteria and data collection modules. If the framework 
can be verified in the future, it could be integrated with 
the other CAD system to optimize the product design and 
work process design. 

In future works, we will apply our results to enhance 
the simulation software such that they will able to produce 
realistic simulations. 
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